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1 ABSTRACT 
O ' 
O 



We report the detection of the linear radical HC4N in the C-rich envelope 
of IRC+10216. After HCCN (Guelin and Cernicharo 1991), HC 4 N is the second 
member of the allenic chain family HC2 n N observed in space. The column density 
of HC 4 N is found to be 1.5 10 12 cm" 2 . The abundance ratio HC 2 N/HC 4 N is 
9, a factor of two larger than the decrement observed for the cyanopolyynes 
(HC2n+iN/HC2n+3N). Linear HC4N has a 3 S electronic ground state and is one 
of the 3 low-energy isomeric forms of this molecule. We have searched for the 
bent and ringed HC4N isomers, but could only derive an upper limit to their 
^ ■ column densities < 3 10 12 cm -2 (at 3-a). 
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1. Introduction 

The detection of long linear carbon-chain molecules in interstellar and circumstellar 
clouds has come as a surprise. Ab initio calculation show that long linear species are usually 
more energetic than their ring or ring-chain isomers and are observed to be less stable in 
the earthbound environment. Yet many linear chains, such as the cyanopolyynes HC2 n +iN, 
are found to be widespread and abundant in the cold circumstellar and interstellar clouds. 
This abundance of long chains and the scarcity of rings in the cold UV-shielded clouds yield 
keys to the formation of large molecules in interstellar space and, may be, can help us to 
understand to the origin of the diffuse interstellar bands (Tulej et al. 1998; Maier, Walker 
and Bohlender 2004). 
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Carbon chain molecules can be classified into several families, depending on the struc- 
ture of their linear backbone (acetylenic backbone with alternating single and triple carbon- 
carbon bonds, or cumulenic backbone with double bonds), on their electronic ground state 
(open shell or closed shell) and on their end groups (H, CN, CH 2 or CH 3 ). The two most 
widespread families are the polyynes (HC 2n+ iN and CH3C 2n+ iH) and the polyacetylenic rad- 
icals (C n H). The most abundant is presumably the family of polyacetylenes (HC 2n H), which, 
unfortunately are non-polar, hence not detectable at radio wavelengths. Note, however, that 
these species have been detected in the mid infrared with the Infrared Space Observatory 
(Cernicharo et al. 2001a,b). 

It has been proposed that the polyynes and carbon-chain radicals form directly in the 
gas phase through reactions of CCH with polyynes, polyacetylenes and/or polyacetylenic 
ions (e.g. HC 2n +iN + CCH -»• HC 2n+3 N + H) (Herbst 1989; Cherchneff, Glassglod and 
Mamon 1993; Millar, Herbst and Bettens 2000). Such reactions tend to insert two triply 
bonded carbon atoms into the carbon backbone, uncoupling the formation of chains with an 
odd number of C-atoms from those with an even number, which could explain the alternance 
between high and low abundances observed as the length of the C backbone increases. For 
example, in the molecular shell around IRC+10216 and in the dark cloud TMC1, the two 
main astronomical sources of carbon chain molecules, the abundance of the C„H radicals with 
an even number of C-atoms is ~ 30 times larger than that of radicals of similar size with 
an odd number of carbon atoms. In contrast, the abundance decrement between successive 
species within the even (or odd) number of C-atoms families, C n H/C n+2 H, is found to be 
only 4-6 (Cernicharo et al. 1987; Guelin and Cernicharo 1991). 

The same seems to apply to the chains terminated by CN: the abundance decrement 
between cyanopolyynes, HC 2n+ iN/HC 2n+3 N, is between 3 and 5 (for n=l-4), whereas HCCN, 
the first member of the 'even' chain family, is found to be ~ 200 times less abundant than 
HC 3 N (Guelin and Cernicharo 1991). So far, this similarity between polyynes and carbon 
chains, could not be further explored as no other member of the HC 2n N was observed in 
space. In this Letter, we report the detection in IRC+10216 of the next member of this 
family, HC 4 N, and show that the HCCN/HC 4 N abundance ratio is 20 times smaller than 
the HC 3 N/HCCN ratio. 



2. Observations 

Contrary to cyanopolyynes which are linear and have x £ electronic ground states, the 
lowest energy form of HCCN is a quasi-linear triplet (Saito et al. 1984, Rice and Schaeffer 
1987). Similarly, ab initio calculations (Ikuta, Tsuboi and Aoki 2000; Aoki and Ikuta 
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1993; Aoki, Ikuta and Murakami 1993) predict that one of the 3 low-energy isomers of 
HC4N is a linear triplet with a 3 E ground state (hereafter HC4N, see Fig. 1), the other 
two being a C3-ring (C3-HC4N) and a bent carbene structure (b-HC4N). The 3 isomers have 
been recently observed in the laboratory (Tang, Sumiyoshi and Endo 1999; McCarthy et 
al. 1999a,b); their microwave spectra are now fully characterized, opening the way to 
unambiguous identification in astronomical sources. 

The astronomical observations were carried out between 1995 and 2003 using the 30-m 
IRAM radiotelescope at Pico Veleta, Sierra Nevada (Spain). Most of the data were taken 
during a sensitive 3-mm wavelength (80-115 GHz) line survey of IRC+ 10216, which has an 
r.m.s. noise of few mK per 1 MHz-wide channel. These data were complemented in August 
2003 with more sensitive observations at the frequencies of SiNC (Guelin et al. 2004). 

We detected a number of U-lines, 9 of which could be grouped into 3 close triplets 
with similar line intensities. The center frequencies of the triplets were harmonically related, 
implying a common linear carrier with a rotation constant B Q = 2302 MHz. The rotation 
constant suggested a molecule similar to C 5 H (B Q = 2287 MHz) and the triplet structure an 
even number of electrons. This pointed out to HC 4 N, and we could quickly check that the 
frequencies of our 9 U-lines agreed with those predicted from the spectroscopic constants of 
Tang, Sumiyoshi and Endo (1999). That HC4N was the carrier of our new lines was definitely 
confirmed when we succeeded in detecting a fourth line triplet (see Fig. 2e and Table 1), 
during a dedicated search made in October 2003. 

The observations were made in the wobbling mode, with the secondary nutating at a 
rate of 0.5 Hz, in order to achieve very flat baselines (see Fig. 2 and below). Only a zero- 
order baseline with an offset of ~ 40 mK, corresponding to the thermal dust emission of the 
circumstellar envelope, has been removed from the spectra. Two SIS 3-mm receivers, with 
orthogonal polarizations and system temperatures of 100 — 130 K, were used in parallel with 
two 1.3-mm receivers. The 3-mm data were co-added, so increasing the effective integration 
time and decreasing the r.m.s. noise. In the case of our best spectrum (Fig. 2a), the 
effective integration time and noise temperature were equal to 52 h and 0.5 mK per 1-MHz 
wide channel, respectively. 

The temperature scale was calibrated with the help of two absorbers, respectively at 
room and liquid nitrogen temperature, using the atmosphere transparency model ATM de- 
veloped by Cernicharo (1985) currently installed in the 30-m telescope on-line calibration 
software (see also Pardo, Cernicharo and Serabyn 2001). Pointing and focus were regularly 
checked on planets and on the strong nearby quasar OJ 287. Between the pointing sessions, 
we monitored the shapes of the strong lines that were observed in parallel with the 1.3-mm 
receivers. The cusped shapes and the line horn-to-center intensity ratios depend indeed crit- 
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ically on the accuracy of the telescope pointing and focus. Thanks to these precautions, the 
intensities of the lines observed several times in the 1995-2003 period were found constant 
within 10%. 



3. Results and Discussion 

Figure 2 shows the spectra covering the 4 triplets, which correspond to 4 successive 
rotational transitions of HC 4 N (from N=18 — > 17 to 21 — > 20). The spectrum centered at 
82.9 GHz, which was observed while searching for SiCN (Guelin et al. 2000), SiNC (Guelin 
et al. 2004) and HCCN, is shown twice in this figure: full scale, i.e. from 82650 MHz to 
83150 MHz (Fig. 2a) and half-scale, around the HC 4 N N=18^17 triplet (Fig. 2b). The 
r.m.s. noise in this 1 MHz resolution spectrum is only 0.5 mK, which makes it one of the 
the most sensitive spectra ever obtained with the 30-m telescope. 

Forty spectral lines are detected in Fig. 2a. Most are weak and partly blended. The 
blended lines are easily resolved, thanks to their characteristic cusps, sharp edges and con- 
stant width in velocity (this shape results from the uniform expansion of the outer circum- 
stellar envelope). We can identify almost 60% of these lines with the help of the molecular 
line catalog maintained by one of us (JC) and described in Cernicharo, Guelin and Kahane 
(2000). The catalog contains the mm transitions of some 1200 different molecular species. 
Of the remaining 40% unidentified lines, at least half could be tentatively assigned to some 
heavy species (see below). The carriers of the unambiguously identified lines are mostly 
12 C and 13 C isotopomers of long carbon chain molecules (C 4 H, C 7 H, C 8 H, H 2 C 3 , H 2 C 4 , and 
HC5N), as well as silicon compounds, such as SiC 4 , SiCN, and SiNC. For the first time, 
we detect a doubly-substituted 13 C isotopomer of cyanoacetylene: H 13 C 13 CCN (see Fig2e, 
96623 MHz; several lines of the three doubly-substituted 13 C isotopomers of HC3N have been 
detected in our 3mm line survey of IRC+10216). 

The number of unidentified lines at the 3 mK level detected in the 4 frequency bands 
is ~15 lines/GHz. Several of these lines probably arise from vibrationally excited C5H and 
CeH and their 13 C isotopomers, species for which we have no accurate laboratory frequencies 
(for example, U82924 and/or U82938 and U82995 could well correspond to the J = 61/2 -> 
59/2 lines of 2 Tl 1/2 CCCCC 13 CH; U82677, U82737, U92034 and U92045 could correspond to 
different £ and n vibronic states of the lowest energy bending mode of C 6 H). Others probably 
come from 13 C and 29 Si isotopomers of of SiC 4 (for example the J=29-28 and J=32-31 lines 
of 29 SiC 4 at 87588 and 96648 MHz respectively and the J=27-26 line of SiC 13 CCC at 82704 
MHz). Finally, some unidentified lines are close to transition frequencies predicted for SiCs 
(U82746, U82876 and U82995), C 5 S (U83038) and KCN (for which we have detected several 
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additional lines that will be published elsewhere), and C 9 H (U83011). Obviously, the U-line 
density is too high at present to claim the detection of any of these species, which all have 
a rich mm spectrum, on the basis of just one, two, or even three weak lines. 

The HC 4 N lines shown in Fig. 2 are cusped (see in particular Fig. 2a, where the lines 
have the highest signal-to-noise ratio). This implies that this molecule is mostly confined 
into the outer part of the circumstellar envelope, presumably in the 15" -radius hollow shell 
where mos t of the free radicals and carbon chains are found (Guelin, Lucas and Cernicharo 
1993). From the intensities of the four triplets we derive a rotational temperature of 25±4 
K, close to that derived for HC3N. 

The permanent dipole moment of HC4N has been calculated by Ikuta, Tsuboi and Aoki 
(2000) to be 4.3 D. Using this value and assuming thermal equilibrium at 25 K, we derive a 
column density for HC4N of 1.5 10 12 cm -2 . 

We have re-calculated the rotation temperature and the abundance of HCCN by adding 
to the 2- mm lines reported by Guelin and Cernicharo (1991) the lines observed in our new 
3-mm survey. We find a rotational temperature of 15±2 K, slightly higher than our previous 
estimate (12±2 K). The HCCN lines are also cusped and the HCCN column density, assuming 
a dipole moment of 3 D (Guelin and Cernicharo 1991), is 1.4 10 13 cm -2 . We note that the 
increase in rotation temperature between HCCN and HC4N is not surprising as a similar 
increase is observed between the cyanopolyynes HC 3 N and HC 5 N. 

The HCCN/HC4N abundance ratio is found equal to 9. This is about twice the decre- 
ment observed between the successive members of the cyanopolyynes family, HC3N/HC5N, 
but 20 - 40 times less than the HC3N/HCCN (~ 200) and HC5N/HC4N (400) abundance 
ratios, which confirms that, like the carbon chain radicals, the even and odd HC n N families 
have parallel, but distinct formation paths. We have also searched in our 3mm line survey 
for the linear triplet HCeN, but without success. The 3cr upper limit to the column density 
of this species, assuming a rotation temperature of 25 K, is 1 10 13 cm -2 , not low enough 
to constrain significantly the HC4N/HC6N ratio, but already 25 times lower than the HC7N 
column density (Cernicharo et al. 1987). 

As noted above, the carbon chain molecules C2 n +iN and HC2 n +iN are thought to be 
formed in IRC+10216 from the reaction of radicals such as CCH and CN with carbon-chain 
molecules or radicals (Cherchneff, Glassglod and Mamon 1993; Millar, Herbst and Bettens 
2000). Similarly, we could consider that HC4N may form through the reaction of C3N and 
C 3 H with HCN and CH 2 , e.g., C 3 H + HCN -> HC 4 N + H, or may be C 3 N + CH 2 ^HC 4 N + 
H. All the reactants are relatively abundant in the outer layers of the circumstellar envelope 
and we expect, by analogy with similar reactions studied in the laboratory, that at least one 
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of these reactions could proceed rapidly. 

Another path proposed for the formation of cyanopolyynes involves ion-molecule reac- 
tions yielding H 3 C„N + or H 2 C n N + , followed by the dissociative recombination of these ions 
(Glassgold, Lucas and Omont 1986; Howe and Millar 1990). The low abundance of ions 
in IRC+10216 makes this path relatively slow (Cherchneff, Glassglod and Mamon 1993). 
Moreover, in the case of HCCN and HC4N, the reaction of HCN with CH^, leading to 
H 2 C2N + + H 2 (or H 3 C2N + + H), does not proceed (Guelin and Cernicharo 1991), making 
the ion-molecule path even slower and less likely than in the case of cyanopolyynes. 

Finally, a last process, formation on dust grains followed by photodesorption when the 
grains reach the outer envelope layers, has been advocated in the case of IRC+10216 (Guelin, 
Lucas and Cernicharo 1993). Such a mechanism may well produce HCCN, HC4N, and the 
cyanopolyynes observed in the outer shell, but would provide no simple explanation to why 
the abundance decrement within each family seems to be so constant. 

The HC4N lines detected during this study belong to the linear triplet isomer. We have 
unsuccessfully searched for the singlet bent HC4N isomer and for the C3-ring cyclic isomer. 
Those latter have a permanent dipole moment almost as large as the linear isomer: 2.96D 
(At„=2.9 D, fi b =0.5 D) and 3.48 D (^=3.15 D, ^=1.49 D), respectively, versus 4.3 D for 
linear HC 4 N (McCarthy et al. 1999a,b; Ikuta, Tsuboi and Aoki 2000). The upper limits 
to the column density of these species are ~ 3 10 12 cm -2 , assuming the same rotational 
temperature than for linear triplet HC 4 N (25 K). 

The most accurate calculations of the HC4N structure to date (level CCSD(T)/cc- 
pCVTZ) predict that the singlet ring isomer of HC 4 N lies 4.2 Kcal/mol below the linear 
triplet isomer (Ikuta, Tsuboi and Aoki 2000). If this result is confirmed, the non-detection 
of the ringed isomer (and of the bent isomer) would further stress the discrimination against 
ringed structures already noted in the case of carbon chain radicals. 

We would like to thank Spanish MCyT for funding support for this work under grants 
AYA2000-1781, PNIE2001-4716 and AYA2003-02785. 
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Table 1: Observed HC4N line parameters 
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Table 1: The number in parenthesis are la uncertainties on the last digit, derived from least 
square fits; they do not include the 5% calibration uncertainty. b denotes a partly blended 
line. The frequencies in col. 2 are calculated from the spectroscopic constants measured by 
Tang, Sumiyoshi and Endo (1999) in the laboratory. 
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Isomer A (linear triplet HC 4 N ( 3 S)) 
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Isomer B (singlet C 3 -ring ( A')) 
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Isomer C (singlet bent) 



Fig. 1. — The three isomers of HC4N considered in this paper. The structures are taken from 
McCarthy et al. (1999b) and Ikuta, Tsuboi and Aoki (2000). Isomer A is linear triplet 
HC 4 N (HC 4 N), isomer B corresponds to the C 3 -ring structure (c 3 -HC 4 N), and isomer C is 
the singlet bent HC 4 N fb-HC 4 N). 
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